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ABSTRACT 

As  a  typical  free  boundary  problem,  a  Stefan  problem  is  studied  from  two 
analytical  and  numerical  points  of  view.  In  the  first  one,  by  changing  the 
dependent  variable  which  stands  for  the  temperature  distribution,  the  Stefan 
problem  is  transformed  into  a  variational  inequality  (V.I.).  It  is  well  known 
that  V.I.  can  be  approximated  by  a  penalized  problem.  The  second  one  is  the 
method  of  the  integrated  penalty  which  gives  a  new  interpretation  of  this 
penalized  problem.  For  these  different  problems,  existence  and  convergence 
theorems  are  given  and,  moreover,  numerical  methods  to  solve  them  are 
presented.  Finally  some  numerical  results  are  given. 
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SIGNIFICANCE  AND  EXPLANATION 

u 

We  can  observe  many  phenomena  involving  Free<  Boundaries  in  various  fields 

V-p  " " 

of  engineering  and  applied  sciences,  for  example,  free  boundary  problems  in 
optimum  design,  the  pollution  of  air  and  water,  the  equilibrium  of  plasma. 

For  such  problems  it  is  important  to  develop  reliable  computational  methods 
which  are  practically  efficient  in  applications.  Naturally,  a  crucial  point 
in  numerical  methods  for  free  boundary  problem  is  how  to  deai  with  the  moving 
boundary.  In  order  to  meet  this  difficulty,  various  approaches  have  been 
proposed;  they  can  be  classified  into  two  groups.  One  of  them  is  to  follow 
the  free  boundary  directly.  The  other  is  to  tranform  the  original  problem 

n ■- 

into  an  auxiliary  problem  with  a  fixed  boundary.  Qu»  technique  uses  the 
penalty  method  and  falls  into  the  second  group. 

In  this  paper  we  apply  this  method  to  the  Stefan  problem  which  arises  in 
the  analysis  of  melting.  of  ice  adjacent  to  a  heated  body  of  water.  The 
essential  idea  of  ouf  method  is  to  transform  the  Stefan  problem  into  an 
initial-boundary  value  problem  for  the  heat  equation  defined  in  a  cylindrical 

domain,  occupied  jointly  by  water  and  ice,  with  an  artificial  heat  absorbtion 

■■n 

*  it:  f 

in  the  ice  region.  Ovf  formulation  is  closely  related  to  the  approach  using 
variational  inequalities. 


\ALL 

The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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NUMERICAL  ANALYSIS  OF  A  STEFAN  PROBLEM 


*  ** 

H.  Kawarada  and  C.  Saguez 


INTRODUCTION 

Free  boundary  problems  appear  in  various  fields  of 
engineering  and  applied  sciences,  for  example,  problems 
in  mechanics  of  continuous  media,  the  equilibrium  of 
plasma,  the  pollution  of  air  and  water  and  others. 

Here  we  restrict  our  interest  to  systems  of  Stefan  type 
where  exists  a  change  of  phase  (solidification,  lique¬ 
faction,  sublimation...) 

So,  we  consider  a  problem  of  Stefan  type  (Pr) Q  in 
one  dimensional  space.  For  this  problem,  we  introduce 
successively  i)  a  variational  inequality  (V.I.)  by 
changing  the  depending  variable  of  (Pr)g?  ii)  a  first 
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penalized  problem  (Pr)^,  which  is  directly  introduced 
from  (Pr)g  by  means  of  the  method  of  integrated  penal¬ 
ty;  iii)  a  second  penalized  problem  (Pr)2  by  changing 
the  depending  variable  of  (Pr)^,  which  arises  as  a 
penalized  problem  associated  with  (V.I.).  The  plan  is 
following: 

Introduction 

1.  Formulation  of  the  problem  and  associated  varia¬ 
tional  inequality 

2.  Introduction  of  the  penalized  problems  (Pr)  ^  and 
(Pr)2 

3.  Some  results  of  existence,  uniqueness  and  conver¬ 
gence 

4.  Proofs  of  Theorems 

5.  Numerical  results 

6.  Conclusion 

1.  FORMULATION  OF  THE  PROBLEM  AND  V.I. 

We  consider  a  one  phase  Stefan  problem  in  one  dimen¬ 
sion.  0 (x, t)  denotes  the  temperature  of  the  solid 
and  x  =  s(t)  is  the  equation  of  the  free  boundary. 

For  simplicity,  we  take  all  physical  constants  equal 
to  1. 

The  problem  (Pr)  n  is  to  find  {9  (x,t)  ,  x  =  s(t)} 
such  that: 

(1.1) 

(1.2) 

(1.3) 

(1.4) 

(1.5) 


3£  329 

3t  ' 

0(0, t)  =  f(t) 

9  (s  (t)  ,  t)  «*  0 
9(x,0)  =  9q (x) 

If  “  _  lx  (s(t)»t)  +h(t)  0  <  t  <  T, 


0  <  x  <  s(t) ,  0  <  t  <  T, 

0  <  t  <  T, 

0  <  t  <  T, 

0  <  x  <  b. 
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(1.6) 


s(0)  =  b. 


where  f(t),  h(t),  0Q(x)  and  b  are  given.  The  term  h(t) 
in  (1.5)  can  represent,  for  example,  a  heat  source 
along  the  free  boundary.  Such  a  phenomenon  can  be  ob¬ 
served  on  the  heat  shield  of  a  space  capsule  during  the 
entry  phase  in  the  atmosphere  or  in  the  continuous 
casting  process  taking  the  convection  in  the  liquid 
into  account  {8]. 

Here  we  assume: 

f(t)  €  C1,t(0,T)*  (0  <  t  <  1) ,  0<f(t)<a 

in  10, T], 

h(t)  €  C1  (0,T)  ,  0<h(t)<8**  in  [0,T]  , 

0  <  b  <  +», 

0o(x)  6  C2(0,b)  ,  0o(b)  *  0,  O<0q(x)<y 

in  [0,b] , 

f  (0)  -  0 0 (0)  . 

Let  Rq  be  such  that  s(t)  <  RQ  for  any  t  £  (0,T) . 

We  denote  by  fr(x,t)  the  extension  of  0(x,t)  by  zero  on 
the  region  (s(t)  <  x  <  Rq,  0  <  t  <  T) .  We  introduce  the 
new  variable  y(x,t): 

rR 

(1.7)  y (x,t)  =  °$U,t)d£. 

'  x 

*  Cm'T(0,T),  m  integer,  0<t<  1,  is  the  space  of  m 
times  continuously  differentiable  functions  in 
(0,T) ,  such  that  the  m  th  derivatives  are  Holder 
continuous  with  exponent  t. 

**  The  condition  h(t)  >  0  can  be  relaxed. 
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Noting  that,  by  the  maximum  principle,  &(x,t)£0,  we 
deduce : 

(1.8)  -  +  h(t)  £  0  in  O^t^T. 

Then,  after  some  calculations,  we  prove  that  y  is 
solution  of  the  following  variational  inequality 
(V.I.):  ((11]  C.  Saguez) 

To  find  y  €  L2 (0,T  ;  H1 (fi) ) , 

6  L2(0,T  ;  L2(n)) , 

S(t)  €  H1 (0,T) , 

fl  =  (0,RQ) 

such  that 

(1.9)  (§£  ,  <D-y)  +  a(y,<)>-y)  >  (-  |f  +  h(t)  ,  $-y) 

+  f  (t)  (4>(0)  -  y  (0,  t) )  , 

4>  6  K  »  {<)>  6  H1  (0)  ;<)»{R)=0,  <j>>0  a.e.  in  ft). 


(1.10) 

y(x,0)  =  yQ(x)  = 

f  °%0(V<K, 

•  x  w 

(1.11) 

s(t)  =  Inf  {x| 

y  (x,t)  =  0}  ,  0  <  t  <  T, 

s (0)  =  b 

(1.12) 

S ( t)  <  RQ 

where  (u,v) 

■  [  °u(x)v(x)dx, 
J  0 

.  .  fR0  du  dv  . 

a(u'v>  *  J0  ara*  dx 

and 

&0(x)  »  |60(x) 

0  <  x  £  b. 

b  <  x  ^  Rq. 

A 


This  inequality  differs  from  an  ordinary  variational 

ds 

inequality  by  the  presence  of  the  term  -  +  h(t)  in 

the  right  member  of  (1.9)/  together  with  s(t)  which 
is  defined  in  (1.11).  An  inequality  of  similar  type 
was  studied  by  A.  Friedman-D.  Xinderlehrer  [3]  and 
Nguyen-Din-Tri  [10],  A  complete  study  of  (1.9) -(1.12) 
is  presented  by  C.  Saguez  (111 . 


2.  INTRODUCTION  OF  THE  PENALIZED  PROBLEMS  (Pr) . 

AND  (Pr)  ^ 

2.1  The  method  of  integrated  penalty 

Here  we  give  the  principle  of  the  method  of  integrated 
penalty  introduced  in  [6]  by  one  of  the  author. 
Originally,  the  penalty  method  has  been  widely  used 
in  optimization  problems  with  constraints.  It  is  also 
a  convenient  tool  for  the  study  of  partial  differential 
equations  (see  for  example  J.L.  Lions  [9] ,  H.  Fujita 
and  N.  Sauer  [4] ) . 

Now  let  us  consider  the  following  penalized 
problem  defined  in  Q^,-  (0,R^)  *  (0,T)  : 

for  any  £  >  0, 


,2.1)  -  |x  (X,t).c  inQT< 

(2.2)  irE  (0,  t)  *  f  (t)  0  <  t  <  T, 

(2.3)  ne(R1,t)  =0  0  <  t  <  T, 

(2.4)  ire  (x,0)  «  ?Q(x)  0  <  x  <  Rjy 
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where  x  =  x(x'fc)  *n  (2.1)  is  characterized  by  the  moving 
boundary  x  =  4;(t)  such  that: 


in  Q<°>  =  {  (x,t) 


0  <  x  <  if/  (t)  ,  0  <  t  <  T} 


1  in  Q 


(1)  ... 
T 


We  assume  that  x  =  ^ ( t )  is  sufficiently  smooth  in  (0,T), 
0  <  ij/(t)  <  (0  <  t  <  T)  and  1 1)  (0)  =  b.  The  problem  (2.1)- 

(2.4)  has  a  solution  tt£ <=  H2,1(QT)*. 

We  easily  verify  that  as  e -*■(),  tt£  converges  to  a 
function  rr°  strongly  in  H1  (Q^,0^ )  and  converges  to  zero 
strongly  in  H^(Q^^)**.  In  fact,  C2,^(Q^^)  is  the 

unique  solution  of  the  following  initial  boundary  value 
problem: 


(2.5) 

3*°  _  32tt° 

3t  3X2 

in  Q‘°> 

(2.6) 

ir(0)  (0,  t)  »  f  (t) 

0  <  t  <  T, 

(2.7) 

1T°(^(t),t)  =  0 

0  <  t  <  T, 

(2.8) 

TT°(X,0)  =  0  q  ( X  ) 

0  <  x  <  b. 

We  define: 

(2.9) 

pe  =  pe(x,t)  =  (|x 

ttg>  (x,t) 

*  Hn‘'n(QT),  m,  n,  integer,  is  the  Sobolev  space: 
{*eL2(QT),  D^6L2(Qt)  (0<u<m), 

D^€L2(Qt)  (0<y<n)}. 


**  H1(Qt)  =  H1'1(Qt) 


4, 
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(2.10)  q£  «  q£(x,t)  *  j  1  P£(£,t)dC  in  Q^. 

Then  we  have: 

Proposition  2 . 1  Let  e  -»  0 

.  0 

(2.11)  p£-  -^-(^(t),t)-6(x-^(t))  infr'(QT), 

(2.12)  q£  -*  "73f(*(t)  *t)  •  (1-X(x,t))  in^lQj). 

The  proof  for  the  elliptic  case  was  given  in  [6J .  If 
we  repeat  almost  the  same  arguments  as  in  [61 ,  we  can 
show  (2.11) - (2 . 12)  . 

The  term  q  is  now  called  the  integrated  penalty. 

£  0 

qe(i|)(t),t)  approximates  the  flux  of  it  (x,t)  through  the 
moving  boundary  x  =  ip  (t)  . 


2.2  The  first  penalized  problem 

H.  kawarada  and  M.  Natori  [5]  transformed  (Pr) q  into 
an  initial  boundary  value  problem  defined  in  a  cylin¬ 
drical  domain  occupied  jointly  by  water  and  ice  for  the 
heat  equation  with  an  artificial  heat  absorption  in  the 
ice  region;  we  have  the  Problem  (Pr)-^; 


To  find: 


f  0£  «  H2,1(Q<r)  n  C1'0^)* 


se  e  C1(0,T) 

.  Qt  *  (0,R0)x(0,T) 


*  cm,n(QT),  m,  n  integer  is  the  space: 


C(QT),  dJ*6C(Qt)  (0<u<m),  D^6CiOT) 

(0<  Y  <  n)). 


-7- 


such  that:  for  any  e  >0 


(2.13) 

ae£  _  32ee  l  _ e 

* 

in  Qt, 

(2.14) 

eE(o,t)  =  f(t) 

0  <  t  <T, 

(2.15) 

isr'V1’  -  0 

0  <  t  <  T, 

(2.16) 

0E  (x,  0)  =  0g(x) 

0  <  X  <  Rq, 

(2.17) 

7^-  |R° (^ xe  e£)  (5/t)d£  +  h(t) 

0  <  t  <  T, 

(2.18) 

Se(0)  =  b  ;  se(t)  <  Rq 

for  any 

t  <£  [0,T]  . 

(2.19) 


Xe  =  Xe(x,t)  = 


0  in  0<x<s  (t),  0  <  t  <  T 


1  in  s  (t)  <  x  <  Rq,  0  <  t  <  T 


The  presence  of  the  penalty  terra  x£0e  in  (2.13)  can 
approximately  replace  (1.3)  if  e  is  sufficiently  small 
Further  we  should  note  that  the  integrated  penalty  of 
(2.17)  approximates  -  —  (s(t),t)  owing  to  proposition 

oX 

2.1. 


As  a  practical  application  of  this  method  in  the 
two  dimensional  case,  we  refer  to  a  study  of  the  be¬ 
haviors  of  frozen  soil  in  a  neighbourhood  of  an 
underground  storage  tank  of  LNG  (Liquid  Natural  Gas) , 
which  is  kept  at  the  temperature  of  -162°C.  A  program 
of  this  method  was  worked  out  and  was  established  for 
practical  uses  by  I.  Yanagisawa  (13J . 


2.3  The  second  penalized  problem 

Let  us  use  the  similar  transformation  as  rn  (1. 


We  define; 


(2.20,  /<*,t)  = 

By  integrating  the  equation  (2.13,  on  <«.*„>' 


obtain: 


(2.21) 


3yl  =  -  i  [R°Xr(5,t)0E(C.t)<»5 


Wow  we  have: 


(2.22) 


e  1  fH0  _  _  i  fR°x 

A  -e  LXe9  e  lx  e  H 


If  s£(t)  <  *<  V 


(2.23, 


=  ~yclx.*> 


If  0  <  x  <  S  (t) 


(2.24) 


*0  ^dd£l«  iv£  (s£(t)  ,t)  . 


E  J  6 


Further  by  (2.17),  we  have: 

dse  _  _  i.  (  °x 
(2.25)  St  ~  e  ]/£ 


0V  ii—  dx  +  h  (t) 


G  3x 


I  y£(sE(t)  ,t)  +h(t) 

e  J 


Therefore,  as  y 
deduce : 


e  is  a  decreasing  function  in  x,  we 
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(2.26)  Ae  =  i  ye  (sG  (t)  ,t)  (1  -  x£)  +  ^X£y£  (x,t> 

°  ^  ye  (sG  (t)  ,  t)  +  ^  X£  (yC  (x,t)  -  ye  (se  (t)  ,t) ) 

-  -h(t)  -  i(yC  -e(^-h(t)))“. 

Thus  we  have  the  following  nonlinear  problem:  (Problem 


SS£L2) 

To  find 


ye€  c2fl(QT), 

Se  €  c2(0,T) 


such  that: 


(2.27) 


3ye  _  a2ye  .  1.  e  .dse  A. 

~£T  e(y  “e(dT-h))  -~--+h 


at 


-  ds 
dt 


in  Q, 


T' 


£ 

(2.28)  ^-(0,t)  =  -f(t) 

(2.29)  yG (Rfl/ t)  =  0 

(2.30)  yG (x,0)  =  yQ(x) 

(2.31)  ^  i  yE(sG(t),t)  +h(t) 

(2.32)  se (0)  =  b. 


0  <  t  <  T, 


0  <  t  <  T, 


0  <  x  <  Rq, 


0  <  t  <  T, 


3.  SOME  RESULTS  OF  EXISTENCE,  UNIQUENESS  AND  CONVERGENCE 
Theorem  3.1  Under  the  assumption  (Al) ,  (Pr) .  has 
a  unique  solution  {6  ,  s  (t)>  which  satisfies: 
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0e  €  H2,1(Qt)  n  C1,0(Qt)  ; 

< 

se  (t)  e  c1{0,T) . 

Theorem  3 . 2  Under  the  assumption  (Al) ,  (Pr) ^  has  a 
unique  solution  {y  ,  sc(t)}  which  satisfies: 

ye  >o,  ye  e  c2/1(qt) ; 
se  e  c2(0,t) . 

Theorem  3 . 3  Under  the  assumption  (Al)  ,  let  e  -»  0 , 
then  the  solution  {ye  ;  s£(t)}  of  (Pr).,  converges  to  a 
solution  {y°  ,  s^(t) }  of  the  (V.I)  such  that: 

i)  ye  -  y°  weakly  in  L2  (0 ,T  ;  H2  (fl ) )  , 

ii)  se  -*•  s°  weakly  in  H1(0,T), 

se  •»  s®  strongly  in  Ct(0,T)*,  (0  <  t  £ 


Theorem  3 . 4  Under  the  assumption  (Al)  ,  let  e  -*•  0, 
then  the  solution  {8  ;  s  (t) }  converges  to  the  unigu 
solution  of  (Pr) Q  such  that: 


i) 


ii) 


0e  *  1° 


weakly  in  H1  (Q^,)  , 

strongly  in  L2 (0) 

for  any  t  6  [0, 

weakly  in 

strongly  ir»  CT(0,T),  (0£T< 


*  CT(0,T)  =  C°'T(0,T) 
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W|H  H  to  M*-* 


is  the  zero-extension  of  6®  into  Q^. 


Theorem  3 . 5  Under  the  assumption  (A.l),  the  (V.I)  has 
a  unique  solution  y°eC3,1(QT)  and  s° (t) =Inf (x;y° (x, t) =0) 
€  C1(0,T) . 

By  means  of  Theorem  3.3,  we  conclude  that  (Pr^ 
interpreted  as  a  penalized  problem  associated  with  the 
(V.I).  Therefore  we  have  a  parallel  structure,  one 
level  of  which  is  a  penalization  with  respect  to  (Par)  ^ 
and  the  other  level  is  with  respect  to  the  (V.I) . 

Thus  we  have  the  general  scheme: 


Method  of  integrated  penalty 


(Pr)0: 

(1 . 1)  -  (1 . 6) 

4 - 

(Pr)  L: 

(2 . 13)  -(2.19) 

1  y (x, t)  =  fR°$(C,t)dS 

I  ye(x,t)  «|R° 

*  '  X 

V.I: 

(Pr) j ; 

(1. 9)  -  (1. 12) 

<i - 

(2.27)- (2. 32) 
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4 .  PROOFS  OF  THEOREMS 


4.1  Proof  of  Theorem  3.1. 


We  put: 


*o  - b  + 


rb 

loV* 


)  dx  +  C.  (b,T)  +  BT, 


and 


C0R0 

Ce  =  “e  +  6' 


CQ  =  sup(ct,Y) 


The  constant  C^(b,T)  is  defined  in  the  appendix  I. 


Define  a  convex  set  in  C[0,T] 


Kj  =  {9  €  C[0,T]  ;  t>  <  «}>  ( t^ )_<<!)  (t2)  <  Rq 


for  0  <  <  t2  <  T  and  <p  (0)  -  b)  . 


First  let  us  study  some  regularity  properties  of  a 
solution  of  the  initial  boundary  value  problem  (2.13)- 
(2.16),  in  which  the  characteristic  function  x£  is 
replaced  by: 

*s 


Xs(x,t) 


f  0  in  0<x<s(t),  0  <  t  <  T, 


[  1  in  s  (t)  <  x  <  R 


0' 


0  <  t  <  T, 


for  some  s(t)  €  K^. 


For  simplicity,  we  denote  this  problem  by  (Pr) 


Proposition  4 . 1 
f  ies: 


(Pr) _  has  a  solution  6  which  satis- 

■3  — - — ■ — - - 


(4.1) 


8  6  H2,1(QT) 
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(4.2) 


Proof 


e  6  c1'°(QT) 


It  is  easy  to  prove  the  existence  of  the 


solution  6  e  H  '  (Q  ) .  Here  we  focus  on  the  proof  of 

*  fill 

(4.2).  Let  us  introduce  a  function  x  e  C  (R  )  for  some 
6  6  (0 , 6 g ]  (<5q>  0)  such  that: 


for  £  6  R 


X  (O  = 


(C  <  0), 

(0  <  £  <  6)  , 
(6  <  U  . 


Obviously,  x  °  (x-s  (t) )  -+  xs  in  L  (QT)  as  6 -*•  0.  Let  0°  be 

the  solution  of  the  initial  boundary  value  problem,  in 

which  x_  °f  (2.13)  in  (Pr)  is  replaced  by  ^(x-sft)). 

5  2  1  as6  aV 

That  is  to  say,  0°€  C  ' X(QT)  satisfies 

15  5  3  X 

-  —  x  (x-s (t) ) -0  in  Qt.  If  we  refer  to  Theorem  4.1  in 

§4  of  chapter  V  in  [7],  we  have: 


(4.3) 


I  3X  ^l'H5  “  \x  (x2't2)  I 

1  ^'5^IX2.“X2^  +  I  ^l”^ 

for  any  (x^t^  and  (x2,t2)6QT, 


(4.4) 


Max  |^-(x,t)  ]  <  C,. 
(x,t)€QT  3X  6 


Here  we  should  note  that  C^,  Cg  and  a  do  not  depend  on 
5£(0,6q].  By  the  theorem  of  Ascoli-Arzela ,  it  follows 
that  there  exists  a  subsequence,  still  denoted  0^, 
such  that: 
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(4.5) 


(4.6) 


ae5  39_ 

3x  3x 


in  C (Qt) , 

in  C  (Qt)  ,  when  <5  -*•  0, 


which  implies  (4.2). 


I 


Proposition  4.2.  The  solution  0=0(x,t)  of  (Pr)  _ 

"  ■  "  _s 

has  the  following  estimates  in  Q^; 


(4.7) 

(4.8) 

(4.9) 


(4.10) 


0  <  0  < 


C0  in  °T 


,t)  Tdx  <  C. 


dx  dt  <  C3, 


2  dx  dt  C4, 


for  any  t  6 (0,T) , 


where  the  constants  CQ,  C2  and  C3  are  independent  of  e , 
but  C4  depends  on  e.  All  the  constants  Cg,  C2,  C3  and 
C4  are  independent  of  s  (6  K^) . 

Multiplying  both  sides  of  (2.13)  by  and  inte¬ 
grating  in  Qt,  we  easily  obtain  (4.8)  and  (4.9).  Here 
we  used  (I. 10)  in  the  appendix  I.  By  (4.8)  and  (4.9), 
we  have  (4.10).  By  using  the  maximum  principle,  we 
easily  prove  (4.7)  | 

Now  we  define  a  transformation  A  such  that: 


(4.11)  s6K^-*-r«As 

(4.12)  r  -  r (t)  -  b  +  i  [tlfR°(x.e)  CC,T)d^]dT 

e  Jo  Jo  8 

+  [  h(x)dt, 

Jo 
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where  0  is  the  solution  of  (Pr)  . 

s 

Let  us  verify  some  properties  of  A. 

Lemma  4.1.  For  any  s  e  K^, 

(4.13)  As 

Proof  First,  we  easily  prove  that  r(t)  satisfies: 
b  <  rlt^)  <  r(t2)  for  0  <  t1  <  t2  <  T 

because  0  >  0  and  b  >  0. 

Next,  by  an  integration  of  the  equation  satisfied  by  6^ 
in  Qfc,  we  obtain: 

(4.14)  (RV(x,t)dx  +  ^  [t([R°{x5e5)  (x,T)dx]dx 

Jo  £  J 0  J o 


Note  that  the  last  term  of  the  right  side  of  (4.14)  is 
estimated  as  follows: 

aft6 

(4.15)  Max  |^-(t,0)|  <  C. (b,T)  for  any  6e(0,5_]. 

°<t<T  dX  1  0 

The  proof  of  (4.15)  is  given  in  the  appendix  I.  Let 
6  0  in  (4.14),  we  have: 

(4.16)  £  Iq^I  °^XS9) (x,t)dx]dT <  |  0q (x)dx + (b,T) 

for  0  <  t  ,<  T. 

Then,  by  (4.12)  and  (4.16),  we  have  (4.13).  | 


Lemma  4.2.  A  is  continuous  in  K^. 

Proof  Let  e (resp.  0')  be  solutions  of  (Pr)„ 

"  •  ■  5 

(resp.  (Pr)  ,)  for  s,s'€K.  and  let  r  =  As  and  r'=»As'. 

S  X 

He  have: 


(4.17) 


!r-r'lctO,T]  i  — -Tls-S1 


C10,T] 


+  rRoie~e,iC(QT)- 


On  the  other  hand,  we  have: 

(4.18)  I«-9'Ic(5t,< 


C0/5R^ 


/l  s~s*  lc(0,Tj  * 


The  proof  of  (4.18)  is  given  in  the  appendix  II.  With 
(4.17)  and  (4.18),  we  easily  obtain  the  continuity  of  A 
in  Kr  | 


Lemma  4.3.  A  is  compact  in  . 

Proof  Prom  (4.12),  we  have: 

(4.19)  |  |^°(XS9)  (s,t)dx  + h(t)  in  0<  t<T. 
Then: 

(4.20)  0  -  a!  -  Ce  in  0<t<T, 

Which  implies  the  compactness  of  A  in  K^.  | 

Proof  of  Theorem  3.1  is  a  bounded,  closed  convex 

in  C(0,T].  By  the  Lemmas  (4.1) -(4. 3),  we  see  that  A  is 
continuous  compact  from  into  itself.  Then  we  can 
use  Schauder's  fixed  point  theorem.  So  there  exists, 
at  least,  one  fixed  point  s(  such  that: 
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(4.21)  s  s  As  in  K^. 

Thus  (Pr)  ^  has,  at  least,  one  solution  which  sat¬ 
isfies  0  6  H2,1(Qt),  0  6C1,O(Qt)  and  x  =  s(t)  6  C1(0,T)  . 
The  uniqueness  of  the  solution  of  (Pr)^  will  be  proved 
at  the  end  of  the  proof  of  Theorem  3.2.  jj| 


4 . 2  Proof  of  Theorem  3.2. 

Let  {0 (x, t) ,s (t) }  be  the  solution  of  (Pr)^.  If  we 
refer  to  the  paragraph  2.3,  then: 

(4.22)  y  =  y (x,t)  *  [R°0(5,t)d£  >  0 

J  x 

and  x  =  s(t)  satisfy  (Pr)~.  Because  of  y  6  C2' 0 (Q_)  and 

1  L 
s  e  C  [0,T] ,  we  see  that: 

(4.23)  §£ec(QT). 

Furthermore,  we  have: 

<4-241  ~f  =  c  l^ls<t,'t)al  +  r  lt<s<t,'t>  +H 

at 

in  0  <  t  <  T 

which  implies: 

(4.25)  s  €  C2(0,T) . 


Therefore  (Pr)2  has,  at  least,  one  classical  solution 
{y€C2,1(QT);  x=  s(t)  6  C2(0,T)}. 

Suppose  that  (Pr)2  has  two  solutions  {y^;si(t)> 
for  i»l,2.  Let  us  estimate  the  difference  between 


ds 


-Z7T  and 


ds 


2 

—  • 

Ja.  f 


we  have : 


(4.26) 


i-jr  — *  c‘lyi*si*t* •  y2(si(t)  >t]  I 

+  |ly2(si  (t) ,t)  -  y2 (s2<t) ,t) | 


♦  i|*i 

e  |  3x 


c(qt) 


_  *  I si  "  s5 1 

C(QT)  x  C(0,T] 


Because 

(4.27) 

we  have 

(4.28) 


ax 


C(Qt) 


^  V 


d3l  ds2 
dt  dt 


<  ^lyi  -  y2l  _ 

C[0,T]  e  C(Qt) 


+  T|sl"s2l 

£  1  z  C  [0  ,T) 

We  have  the  following  inequality  (see  the  appendix  m) : 

rT 


(4.29) 


+  2 


T 

ds^  ds2| 
•HT  ’  dTj 


C[0,T] 


)dt. 


Prom  (4.29),  we  obtain: 
(4.30) 


2TQe  ids.  ds2 

Iyi  -  y2l  -  -  “  rHF  "  "dt" 

1  C(Q_  )  E  T0  1 
T0 


C[0,T01 


for  0  <  TQ  <  c . 
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Using  (4.28)  together  with  (4.30),  we  have 


(4.31) 


ds^  ds2 

trt  "  ■ar 


cio,tq] 


.  ,  2T0  .  C0T0.  dsl  ds2 

i  «i=fj  — 1  -  -at-  c|#  , 


for  0  <  Tq  <  e . 


Choose  Tq  *=  Tq  (e )  to  satisfy: 


(4.32) 


That  is: 


2T0  C0T0 

0  <  —  -  -  +  — - —  <  1  for  any  e  >  0. 

o  e 


(4.33)  0  <  TQ(e)  <  Inf  +  CQ  - /cj  +  2CQ  +  9)  ,e  }  . 

Then  we  have  a  contradiction  and  we  deduce 
ds.  ds, 

(4.34)  -3F(t)  =  inO<t<TQ(e). 

Because  s^tO)  =  s2(0)  =  b,  we  have 


(4.35) 


s1(t)  =  s2(t) 


in  0  <  t <  TQ(e) . 


With  (4.33)  and  the  same  procedure  {[ 
We  obtain: 


f^rrr1  +11  tines' 


(4.36) 


S^t)  =  S2(t) 


in  0  <  t  <  T , 


which  implies: 


(4.37)  y^Xjt)  -  y2(x,t)  in  Q^. 

Thus  this  proves  the  uniqueness  of  the  solution  of 
(Pr)2*  Prom  this  uniqueness,  we  easily  obtain  the 
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uniqueness  of  the  solution  of  {Pr)^.  Let  (8^  ,  s^(t)} 
and  {8 j  »  Sj(t)}  be  two  solutions  of  (Pr)^.  Then 


*1 


/t)d£ 


and 


*2 


U,t)d£ 


satisfy  (Pr)_.  From  the  uniqueness  of  the  solution  of 

£  V 

(Prjj#  we  deduce: 

(4.38)  s1(t)  =  s2(t)  in  0,<  t^T. 

On  the  other  hand,  we  have, 


(4.39) 


($/t)  —  0  2 


(5,t)d£  =  0 


for  any 


(x,t)  e  Qn 


If  we  differentiate  (4.39)  with  respect  to  x,  we  have 
(4.40)  O^Xjt)  «  0  2  (x ,  t )  in  Qt. 


4.3  Proof  of  Theorem 

Here  we  denote  by  {ye  ;  se(t)}  the  solution  of 
(Pr) 2* 


Lemma  4.4. 

yE  verifies 

the  following  estimates: 

(4.41) 

-cft  <  <  o 

0  *»  3x  ~ 

and  0  <  ye  <  CQR0  in  QT 

(4.42) 

2dx  <  C-  for  any  t  6(0, T) 

'0  1  3x* 

“  a 

(4.43)  ||  |^~|2dxdt<C 
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r 


(4.44) 


II 


IX 


0T'  3x 


dx  dt  <  C4 . 


We  easily  obtain  (4.41)-(4.44)  from 

(4 . 7)  (4 . 10)  in 

Proposition 

4.2. 

Lemma  4.5. 

se (t)  verifies  the  estimate: 

(4.45) 

|dse  <  c 

nst*  2  i  C7 ' 

»  LZ10,T)  1 

where  C,  is 

independent  of  e . 

1 

Proof  From  (2.27) ,  we  have: 

(4.46) 

=  h.3zi  +  sV 
^  3t  ^ 

in  0  <  x  <  se  (t)  , 

0  <  t  <  T. 


From  (4.46),  we  have 


3 1  i 


dx  dt 


|f 

Qrp  3* 


dx  dt' 


+  R0-T-Bz]. 

By  (4.42),  (4.43)  and  (4.47),  we  deduce: 

.£ 


(4.48) 


Lemma  4.6. 
5  0 


If  e  -»  0,  then  {yE ,x  ■  sG  (t)  )  converges  to 


{y  »s  ( t) }  in  the  following  topology: 
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(4.49) 

E  0 

s  -*-  s 

2 

strongly  in  L  (0,T) , 

strongly  in  Ct(0,T),  (O^t^j) 

(4.50) 

dse  ds° 

dt  dt 

weakly  in  L^(0,T), 

(4.51) 

ye  -*>  y°(>o) 

2 

strongly  in  L  (QT) , 

(4.52) 

izi  h. 

3x  3x 

2 

strongly  in  L  (QT) , 

(4.53) 

ax^  ax 

weakly  in  ( 0 , T  ;  ( 0 , Rq ) ) , 

(4.54) 

izl  ^  iZ 

at  at 

weakly  in  (0,T  ;  (0,Rq) ) • 

Define  a  convex  set  K„,  such  that 

K  =  { $  e  L2  (fi )  ;  <{>  >  0  a .  e .  in  Q }  . 

If  we  multiply  (2.27)  by  (4>-ye),  $eK2,  and  integrate 
in  Qt»  then  we  have: 

(4.55)  [  (^- ,  4»-yE)dt  -  [  (^“V  »  $-ye)dt 

J0  3t  Jo  ax'6 

-T  ,T 

-if  ((ye+ea  )”  ,  ij>~ye)dt  «  (a  ,<J>-y  )dt 

eJ  o  e  ;0 

€  K2 . 

dse 

Where  a^  =  -  -g^-  +  h(t)  . 

Now,  using  the  monotone  property  of  (ye  +  eaE)  , 
and  because  $6K 2»  we  have: 


(4.56) 


-  (ye+ea£)  ~  ,<p  -  y£) 

=  £■(♦”-  (ye  +e )  ■" , <t>  -  (ye+cae) ) 

-  ( (ye+ea£)“,a£) 

<  -( (y£+ea£) ",ag) . 

Using  (4.45),  we  have 

2 

(4.57)  ea£  -*•  0  strongly  in  L  (0,T)  as  e  -*■  0 

By  (4.51)  and  (4.56),  we  have 

(4.58)  (ye+ea£)”  (y°)“  =  0  strongly  in  L2  (QT) 

as  e  -*■  0 

Prom  (4.58),  we  have 

(4.59)  f  ((ye+ea  )~,a  )dt  -*■  0  as  c  ■*  0  . 

J  0  EE 

Therefore,  using  Lemma  4.6,  (4.56)  and  (4.59),  we 
obtain  at  the  limit  in  (4.55): 

fTa°  n  rT  .2  0  n 

(4.60)  (-^r-  >  4>-yU)  dt  -  (£_X_  ,  4>~yU)  dt 

J0  c  JO  dx 

rT 

>  j  (aQ,<j>-y0)  dt. 

o 

Here  +  h.  The  partial  integration  of  the 

second  term  of  the  left  side  of  (4.60)  yields  (1.9). 
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4.4  Proof  of  Theorem  3.4. 

For  simplicity,  we  assume  f  =  0  in  (0,T) . 

Step  1;  convergence  of  0£  and  sc. 

We  already  proved  the  convergence  of  s  (t)  in 
Lenuna  4.6.  From  Proposition  4.2,  we  can  extract  a 
subsequence,  still  denoted  0£,  such  that: 

(4.61)  0e  --  w  weakly  in  H1  (QT>  , 

9e  -*■  w  strongly  in 

and  we  have  ((4.8)  Proposition  4.2): 

(4.62)  f*°l|5<--t)|2dx  <  C2  Vt£(0.T). 

Define: 

fl°(t)  *  (x  ;  0 <  x < s°(t) }  (0  <  t  <  T) , 

n^t)  -  n  -  n°(t)  , 

q£  -  u  n°(t) , 

1  0< t<T 

and  Qi  -  U  n1(t). 

1  0<t<T 

Lemma  4.7.  The  function  w  satisfies: 

(4.63)  w|  o  6  L2(0,T  ;  hJ («° (t) > ; 

°T 

(4.64)  w I  ,  -  0. 

Iqt 

Proof  First  we  prove  (4.64).  We  have. 


(4.65) 


f|  Xe|6e|2  <*xdt  <  e*C8, 


T 


where  Cg 
we  get 

is  independent  of 

(4.66) 

||  i!w|2dxdt 

qt 

By  using 

(4.62),  we  have 

(4.67) 

|  w(x1,t)  -  w(x 

1 

<  /C^]x1-x2 | 

which  implies,  with  (4.64) 

(4.68) 

w|  o  =  0 

lx=sU(t) 

Since  we  H1(qt),  we  have  (4. 


.  From  (4.59)  and  (4.61), 

0. 

t)|  <  |^||2(..t)|dx 

vx1,x26(0,r0)  ,  vte(0,T), 

Vt  6(0, T)  . 

63)  .  | 


Lemma  4.8. 


(4.69) 


w  g  satisfies : 


jlQ0  5T*d*at  +  jI00H!; 


dx  dt  =  0 


V<t»  6  L2(0,T  ;  Hq  (fi° ( t) ) ) 


Proof  For  any  4> ,  we  can  define  <J>  such  that: 

ii)  $e(0,t)  *  0  (t 6(0 ,T) } ; 


^  ^  9  0  /  /v  %  J  _•  » 


i)  <T  €  H  9  (QT)  ; 


iv)  $e  -*  $  strongly  in  L2  (QT)  and 

-*•  strongly  in  L2  (QT>  when  e -*  0. 

»v>  1 

Here  $  denotes  the  zero-extension  of  <{>  into  QT.  If  we 

multiply  (2.13)  by  <)>e  and  if  v»e  intearate  in  QT,  we  have 

l4'70)  IIqt  tt  ******  +  ||0tC  ^rdltdt  - 

At  the  limit,  we  obtain, 

(4.71)  f  J  Qj£  *  d* dt  +  idxdt  *  °-  ■ 

Now  let  us  consider  the  following  initial  boundary 
value  problem  (Pr)  Q; 


36  a20 

3t " 

in  qJ, 

6(0, t)  =  0 

0  <  t  <  T 

0 (s° ( t) , t)  *  0 

0  <  t  <  T 

0 (x,0)  =  0Q(x) 

0  <  x  <  b 

1 

0  7 

Here  we  should  note  that  sec  (0,T).  Refermg  to 

Friedman  [2J,  we  easily  obtain  that  (Pr)  0  has  the 

0  2  10  s 
unique  solution  0  £  C  '  (Q_)  which  satisfies 

0  i 

%-<s°(t>  ,t)  €  C  ( 0  ,T)  . 

Lemma  4.9.  w  satisfies; 

(4.72)  w(x, t)  =  9°(x,t)  a.e.  in  qJ. 
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Proof  Put  W  =  w-e.  Then  W  verifies: 

<4-731  |L>  !?*dxdt  ♦  1L>  I2!£ 


dxdt  =  0 


V$  e  L2(0,T  ;  hJ  («°  (t) ) )  . 


Put  $  =  W  in  (4.73).  Then 


(4.74) 


j|W(*,T) |2  + 


IMS 


dx  dt  =  0, 


which  implies  (4.72).  f 

We  modify  w  on  a  null  set  so  that  (4.72)  holds  for 
every  (x,t)  in  Q°,  i.e.,  w(x,t)  =  e°(x,t). 

Step  2;  9°(x,t)  and  s®(t)  satisfy  the  Stefan  condition 

(1.5) .  Define 

(4.75)  qe  =  qG(X,t)  =  j  ° (|  X£  e£ )  U , t) d£ . 


Lemma  4.10. 


(4.76) 


If  e  -►  0,  we  have 
qe-*  -e°(s°(t)/t)  <1-X°(x,t)) 


weakly  in  L  (Q^) . 


Proof  From  (2.12)  (Proposition  2.1),  we  obtain 
(4.77)  qE  -  -9°(s°(t) ,t)  (1  -  X°(x,t)  in  &  (Q^) . 

From  (4.7)  (Proposition  4.1),  it  follows  that: 


(4.78) 


-  h  =  q£(0,t)  >  q£(x,t)  >0  in  Q  , 
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then  we  have: 


(4.79)  |qe|L2(  }  <  /R^.  (C7  +  g/^)  . 

With  (4.79),  (4.77)  implies  (4.76).  | 

Finally  we  show  that  0^  and  s^(t)  satisfy  the 
Stefan  condition  (1.5).  From  (2.17),  we  have: 

(4.80)  ^r{1"Xe>  “  q£U-Xe>  +h(l-xe)  m  qt. 

Integrating  (4.80)  in  (0 <  t <  T)  and  passing  to  the 
limit,  we  have: 

JLo  TRT  ^  |  0  ^(S°(T),T)dxdt 

Qt  Qt 

+  UQohdxdT- 

t* 

Here  we  used  (4.76).  By  Nikodym's  theorem,  we  have 
(4.82)  s°(t)  =  b- ~(s°(x)  ,x)dt  +  |\(T)dT. 

Since  ~-(su(t)  ,t)  6  C[0,T]  ,  (4.82)  implies  (1.5). 
Because  of  the  unique  existence  of  the  solution  of 
(Pr)Q,  all  the  sequence  { 0e (x, t) ,se (t) }  converges  to 
{Q°(x,t) ,s°(t) } . 

4.5  Proof  of  Theorem  3.5. 

From  Theorems  3 . 3  and  3.4,  we  have : 

2 

strongly  in  L  (Q?) 

2 

strongly  in  L  (Q^) . 


3X  3x 


and 


Since 


a.e  in  Q^, 

a.e  in  Q?. 

1 

=  y°(R0,*)  6  h2(o,t). 

Prom  (4.84)  and  (4.85),  we  have 

(4.86)  y°(x,t)  =  fR°9'°(C,t)dC. 

•*  x 

Since  0°  verifies  0<  0°<  CQ  in  Q®,  we  have 

(4.87)  s°(t)  -  Infix  ;  y°(x,t)  «  0}  (0 <  t <  T) . 

Here  x=s°(t)  is  the  free  boundary  of  (Pr)  g  correspond¬ 
ing  to  9°.  {6°  ,  s°(t)}  is  the  unique  solution  of 

(Pr)Q.  Therefore  the  (V.I)  has  a  unique  solution. 

5.  NUMERICAL  RESULTS 

Here  we  present  the  numerical  methods  used  to  solve 
Problem  (Pr) 2  and  the  (V.I).  The  results  concerning 
(Pr)^  were  presented  in  H.  Kawarada-M.  Natori  (5). 

We  use  the  following  algorithms. 

5.1  Problem  (Pr) 2 
1)  Initialization. 

Put  Yq  -  yQ(x),  Sq  «  b  for  n  =  0. 


(4.83) 
we  have 

(4.84) 

We  obtain 


i*l 

dX 


=  -0' 


=  - 


dX 


(4.85) 


0  s  y£ (Rq#  * ) 
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4 


£  X  • 

3)  Determine  y  ^  solution 


2  £  ,  i 

*  r  • 


(5.1) 


yn+i,i-yn  _ a  yn+i .  if..Cfi 
—  T~  c  vyn+l 


At 


dx 


se  ' 1  -  sc 
,sn+l  sn 

At 


-  h 


n+1)r 


e  t  i  e 

sn+l  “  sn  ,  un+l 
— it — th  • 


(5.2) 

(5.3) 


^±1(0)  =  -fn+1. 


dx 


yi;i«Ro>  -  0 


4)  Determine  s^J  by  the  expression: 

, i+1 _  e 
n+1  n  _  1  e  h 

- 1 - e  yn+l{sn+l'  +  h 


n+1 


s)  <;i+1= 

6)  Test  of  convergence; 


(0  <  o>  <  1) 


e  ,i+l 


=  ve'f 


If  the  test  is  verified,  sn+1  *  sn+1  ,  yn+1syn+i 


If  not,  i  =  i+1  go  to  3. 

7)  If  n+1  =  N  (number  of  step  in  t) ,  stop. 
If  not,  n  =  n  +  l  go  to  2. 


The  nonlinear  equation  (5.1)-(5.3)  is  approximated  by  a 
finite  difference  method  (or  a  finite  element  method) . 
This  problem  is  solved  by  an  over-relaxation  method. 


The  exact  solution  is: 


6 (x,t)  -  - 2x  +  t  +  ^)+;  s(t)  =  2  - /3-2t  . 

To  solve  (Pr)2/  we  took  RQ  =  T=1.  We  divided  (0,RQ) 
into  30  steps  and  (0,T)  into  20  steps.  The  results, 
presented  here,  were  obtained  for  e  =  10  ^  and  m  =  0.15. 

The  computer  time  was  approximately  equal  to  4s  on 
IBM  370-168.  For  the  V. I ,  we  took  30  steps  of  space 
and  10  steps  of  time.  The  coefficient  u  was  chosen  to 
be  equal  to  0.3.  The  computer  time  was  the  same  or¬ 
der  as  for  the  preceding  method.  The  results  on  the 
free  boundary  were  very  comparable,  which  are  shown  in  the 
following  table: 


t 

Exact  sol. 

(Pr)  2 

(V.I) 

0 

0.2679 

0.2679 

0.2668 

0.1 

0.3267 

0.3303 

0.3235 

0.2 

0.3875 

0.3956 

0.3845 

0.3 

0.4508 

0.4576 

0.4468 

0.4 

0.5168 

0.5290 

0.5135 

0.5 

0.5858 

0.5977 

0.5832 

0.6 

0.6584 

0.6711 

0.6578 

0.7 

0.7351 

0.7488 

0.7363 

0.8 

0.8168 

0.8325 

0.8196 

0.9 

0.9046 

0.9208 

0.9112 

1 

1.000 

1.023 

0.9965 

We  conclude  a  very  good  agreement  of  the  numerical 
results  with  the  exact  solution. 

On  the  other  hand,  the  method  of  fixed  point  will 
be  used.  Particulary  we  notice  the  utilization  of  the 
method  presented  in  M.  Sermange  [12] , 
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6.  CONCLUSION 

The  methods,  presented  here,  have  a  fairly  general 
character  for  the  free  boundary  problems  appearing  in 
the  system  including  change  of  phase.  They  permit,  in 
particular,  to  take  into  account  heat  source  along  the 
free  surface  and  give  efficient  numerical  algorithms. 
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APPENDIX  I 

1st  step  We  shall  consider  the  following  initial 
boundary  value  problems  (A) ,  (B)  and  (C) : 


( 


(A) 


(B) 


(C) 


39. 

at 


a2e 


3  x 

8  (0,t)  =  f  (t) 

0  (b, t)  -  0 

{  e(x,o)  =  eQ(x> 


in  0  <  x  <  b,  0  <  t  <  T, 

0  <  t  <  T, 

0  <  t  <  T, 

0  <  x  <  b. 


0(0, t)  =  f(t) 


0 (R,t)  =  0 
0(x,O)  =  &Q(x) 


38 

at 


a20 


in  0  <  x  <  R, 
0  <  t  <  T, 
0  <  t  <  T, 
0  <  t  <  T, 
0  <  x  <  R, 


in  0  <  x  <  R,  0  <  t  <  T, 


3x 

9(0, t)  =  f(t) 


9x(R,t)  *  0 


0 (x,0)  =  0Q(X) 


0  <  t  <  T, 
0  <  t  <  T, 

0  <  x  <  R. 


Here  R  in  (B)  and  (C)  is  an  arbitrary  number  which 

g 

satisfies  R>b>0.  s(t)  in  x  (B)  satisfies 

i)  s«C°(0,T)  ; 

ii)  s(t)  is  monotone  increasing  in  t  ; 

iii)  s(0)  »b,  S(T)  <  R. 
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Let  0^,  and  @2  be  the  solutions  respectively  of 
(A),  (B)  and  (C)  under  the  assumption  (Al) .  Then,  from 
the  maximum  principle,  we  have: 


(1.1) 


01(x,t)  <  0  (x,t)  <  02(x,t) 


in  0  <  x  <  b, 
0  <  t  <  T. 


From  (1.1)  it  follows: 

L5 

3x 


(1.2) 


I  Tv"*0' I  0fn 

C  (0,T] 


30. 


<  Sup{  |  -j^(0,  t) 


30. 


I  o  '  I  "a v  (0' t)  I  o  ^  * 

C°tO,T]  dX  C°[0,T] 


30. 


2nd  step  We  shall  prove  that  |  C 0 ,  t) 


uniformly  bounded  for  a:iy  R  (>b)  . 
We  introdu* 

(0  <  t  <  T)  : 


Cu[0,T] 
30 


is 


We  introduce  Green's  function  which  satisfies  g^-(R,t)=0 


2 

N(x,t;£-R,T)  =  —  1  -  exp(  -  -  ) 


2/ir  (t-x) 


.  - 1 -  exp{  -  ■  ) 

4(t-T) 


*  K(x,t  ;  C-R,t)  +  K(x,t  ;  -£+R,t) . 
Integrating  Green's  identity: 

JL  (n  ill  _  e  M )  -  — (n0  )  =  o 
v  35  2  35  3t  21 


on  the  domain  0  <  5  <  R»  0  <  5  <  t  <  t-e  and  letting  e  -*■  0, 
we  obtain: 


rb 

(1.3)  e2(R-x,t)  =  j  N(x,t  ;  £-R,0)90(£)d£ 

ft  302 

-  j^N(x,t  ;  -R,t)  (0,t ) <3t 

+  jt  §f(x,t  ;  -R,2)f  (T)dx. 


Let  denote  v(t)  (0,t).  We  differentiate  (1.3)  with 

respect  to  x  and  x-*-R-0.  We  obtain: 


(1.4) 


(t)  «  2 |  V(T)f£(R,t;-R,T)dT 

+  2 [  G(R,t  ;  -R,t)  f  (T)dT 
JO 

fb 

-  2 j  G(x,t  ;  ?-R,0)60(C)dC, 


where  G(x,t;£,T)=K(x,t;{i,T)-K(x,t;-S,T)  and  f  (t  ?  . 

We  note  that  the  following  inequality  holds: 

.  t  .00 

U.5)  |  (R,t  ;  -R,t)  |dt  <  Aft  ?~e“Ydy 

JO  3X  °J  4R^  /y 

t 

»  Y(R)  , 

where  AQ  is  an  absolute  constant.  Taking  R  sufficient¬ 
ly  large  so  that: 


(1.6) 


we  have 


0  <  y(R)  <  3  * 
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Let  R  -*■+<» .  Then  D(R)  -'•D(+°°)  <  +®.  We  see  that  for  any 
R  which  satisfies  (1.6),  there  exists  an  absolute 
constant  Dq  such  that: 

(1.8)  D (R)  <  Dq. 

Thus  there  exists  some  constant  such  that: 

(1-9)  ! VIC (0/T)  =  D1  f°r  any  R-b* 

30. 

Letting  C.  (b,T)  =Sup{D  ,  |-^(0,t)|  }, 

1  x  ox  C  [0,T] 

we  obtain 

a  a  ^ 

(I. 10)  4£-(0,t)  <  C. (b,T) 

dX  C [0 ,  T]  1 

for  any  R>  b  and  any  5  €  ( 0 , 5  Q  ]  ( 6  >  0 )  . 
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APPENDIX  II 


Putting  0  -  0 '  »  w#  we  have 

,2 


(II. 1) 


|?  -  -  rxsw'F0s*  •{xs"xs')- 

**  ax 

3  w 


Multiplying  both  sides  of  (II. 1)  by  n  and  integrating 


in  Qt#  we  have: 


(II. 2) 


ft[R°!|?1 

J  0-*  o  3t 

C^ftfR 


2  dx  dt  +  ^ 


°||£(x,t>|2dx 
«  dX 


c0fc[  0,  _  I  il2£| 

-joJQ  !xs  xs.llatl 


dx  dt. 


Here  we  used  the  monotone  property  of  s(t).  By  use  of 
Schwartz's  inequality,  (II. 2)  becomes 


(II. 3) 


f fc fR0 |3wi2  ,  ,  +  ifR0|3w 

]Jo  'it'  dxdT  *Jo  '3* 


(x,t) |  dx 


<  -^/t|s  -  s’  |  2(  ) 

J  C[0,T]  L  (Qt) 


for  any  t  e  (0,T) 


From  (II. 3)  we  have 
(II. 4) 

By  (II. 3)  and  (II. 4),  we  have: 


C  _  1/2 

<  -9./ris  -  s'  l 

'3t  t2(Oj)  "  E  C‘°'T) 


(II. 5)  fR0|3i(x,t)|Jd*  <  2(-^)2T|s-s'lcl()/Tr 


On  the  other  hand,  we  have 


(II. 6) 


|w(x,t>|  < 


<  /2RqT*-^-|s  -  s' 


(II. 6)  implies  (4.18). 


-40- 


APPENDIX  HI 


Putting  y ^  -  y2  =  W,  we  **ave 

'  9W  9  2W  1  . 

»  ’  7?  *  ^ 

1+ea1)  "  -  ^(y2+ea2)  +  a1  -  a2  in  QT> 

< 

|^(0,t)  -  0 

0  <  t  <  T, 

W(BQ,t)  =  0 

0  <  t  <  T, 

W(x,0)  =  0 

0  <  x  <  Rq. 

ds . 

Here  a±  =  -  +  h 

(i=l, 2) .  Let  U(x,£  ;  t  ;  t)  be  the 

Green’s  function  of 

the  linear  system; 

'  3u  32u 

3t  ‘  17 

in  Qt, 

1 

|H(0.t)  -  0 

0  <  t  <  T, 

U (Rq» t)  =  0 

0  <  t  <  T. 

Then  we  have 

W(x,t)  =  f  dT 

J  0 

J^°  U (x, t  ;  Z,t)  (|-(y1+ea1)“ 

-  |(y2+ea2)-  +  a.^  -  a2)dC. 

By 

(4.29),  we  have 

(m 

.1)  |W(X,t) 

1 

*  CH 

,R 

0  u(x,t;(,T) {|iyi-y2i +  2 1 ai"a2 ^ 
o 

Noting  that 
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we  have 
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